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In superconducting ferromagnets for which the Curie temperature T™, exceeds the superconduct- 
ing transition temperature Tc, it was suggested that ferromagnetic spin fluctuations could lead to 
superconductivity with p-wave spin triplet Cooper pairing. Using the Stoner model of itinerant fer- 
romagnetism, we study the feedback effect of the p-wave superconductivity on the ferromagnetism. 
Below Tc, the ferromagnetism is enhanced by the p-wave superconductivity. At zero temperature, 
the critical Stoner value for itinerant ferromagnetism is reduced by the strength of the p-wave pair- 
ing potential, and the magnetization increases correspondingly. More important, our results suggest 
that once Stoner ferromagnetism is established, is unlikely to ever be below Tc- For strong and 
weak ferromagnetism, three and two peaks in the temperature dependence of the specific heat are 
respectively predicted, the upper peak in the latter case corresponding to a first-order transition. 



PACS numbers: 71.10.-w, 71.27.-|-a, 75.10.LP 

Due to the strong interplay between conventional su- 
perconducting (SC) and ferromagnetic (FM) states, the 
exploration of their possible coexistence in the same crys- 
tal might have seemed fruitless, but has nevertheless at- 
tracted a great deal of interest recently. This possible 
coexistence was first proposed by Ginzburg more than 50 
years ago Several years later, Larkin and OvchinnikoA*^ 
and Fulde and Ferrell^ independently developed a mi- 
croscopic theory of this coexistence in the presence of 
a strong magnetic field, based upon a spatially inho- 
mogcneous SC order parameter, presently referred to as 
the FFLO state. Meanwhile, Berk and Schrieffer sug- 
gested that conventional s-wave superconductivity in the 
paramagnetic phase above the Curie temperature Tm is 
suppressed by critical ferromagnetic fluctuations near to 
Tm— . However, more recent calculations showed that 
conventional s-wave superconductivity can form in the 
weakly FM regime close to a quantum phase transition^. 
In addition. Fay and Appel predicted that p-wave su- 
perconductivity could arise in itinerant ferromagnets^. 
Their pioneering work indicated that longitudinal fer- 
romagnetic spin fluctuations could result in a p-wave 
"equal-spin-pairing" SC state within the FM phase. 

Experimentally, a major development occurred with 
the observation by Saxena et al. that UGe2, nominally 
an itinerant FM compound, undergoes an SC transition 
at low Tc values under high pressure^.. An SC state was 
also found in other itinerant ferromagnets such as ZrZn2 
and URhGe^i^. In each case, the regime of the SC phase 
appears completely within that of the FM phase, suggest- 
ing a cooperative effect between the SC and FM states. 

These experimental achievements have stimulated re- 
newed theoretical interest in the subject. Recently, a 
large effort has been devoted to the understanding of the 
underlying physics of the coexisting SC and FM states, 
with a focus upon the SC pairing mechanism and the 



orbital symmetry of the SC order parameter. Although 
earlier works by Suhl and Abrikosov suggested that an 
s-wave pairing interaction between conduction electrons 
could be mediated by ferromagnetically-ordered local- 
ized spins, such as by impuritie a^*^'^^ , recent studies of 
these SC ferromagnets^iSii^ have assumed that the itin- 
erant electrons involved in both the FM and SC states 
are within the same ban d-'^^'^'^i^^'^^i"'^^'^^ . Some of these 
studies assumed conventional s-wave pairing. For exam- 
ple, Karchev et al. studied an itinerant electron model 
in which the same electrons are responsible for both the 
FM and SC states^. In that study, the Cooper pairs 
were assumed to be in a spin-singlet state, and the ferro- 
magnetism was described within the Stoner model. How- 
ever, the resulting SC ferromagnetic state was shown to 
be energetically unfavorable when compared to the con- 
ventional, nonmagnetic SC stat^l^. A possible exception 
to this incompatibility could occur if the magnetic in- 
stability were to arise from a dynamic spin exchange in- 
teraction, as discussed by Cuoco et alM. On the other 
hand, a number of other workers avoided the likely in- 
compatibility of the SC and FM states by assuming a 
spin-triplet SC order parameter with p-wave orbital sym- 
metry, for simplicityi^ii^iii. Kirkpatrick et al. indicated 
that a p-wave SC state meditated by ferromagnetic spin 
fluctuations is more likely to coexist within the Heisen- 
berg FM phase regime than within the paramagnetic 
phase regimei^. Machida and Ohmi studied the prop- 
erties of a p-wave SC ferromagnet phenomenologically^. 
More recently, a microscopic model of the coexistence of 
a nonunitary spin-triplet SC state with a weakly itin- 
erant FM state was developed by Nevidomskyy''^^. The 
present nature of the SC coexistent with the FM state 
in these ferromagnetic superconductors is still somewhat 
controversial, although increasingly, additional experi- 
ments on the U-based materials have provided increasing 
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support for a spin-triplet state rather than a spin-singlet 
onei^^i^o^. 

Most theoretical studies have focused primarily on the 
effect of the established ferromagnetism upon the na- 
ture of the coexistent superconductivity, as summarized 
above. However, to fully understand the interplay be- 
tween the SC and FM states when they coexist, one 
should also study the feedback effect of the superconduc- 
tivity upon the ferromagnetism itself, as has been done 
in only one study to date-^. 

Here we study explicitly the effects of the p-wave pair- 
ing on the FM ordering, using the Stoner model of itin- 
erant ferromagnetism as the starting point. We cal- 
culate the critical Stoner parameter Uc, the magneti- 
zation m, and the two parallel-spin p-wave gap func- 
tion magnitudes, A±, respectively, as functions of the 
pair-interaction strength V. We also discuss finite- 
temperature properties, including the T-dependencies of 
these order parameters and the specific heat C(T). 

We take the Hamiltonian for the ferromagnetic super- 
conductor to have the form 



H 



FM+SC 



k,k' 



where cr = ± represent the single-particle spin states, 
and the single-quasiparticle part of H comprises the 
Stoner model for itinerant electrons, where et is the non- 
magnetic part of the quasiparticle dispersion, fi is the 
chemical potential and M = U{{n+) — (n_))/2 is the 
magnetic molecular-field with U the Stoner exchange in- 
teraction, and V is the sample volume. The pairing po- 
tential is taken to have the p-wave form—,, Vsc(k, k') = 
— Vh • k'. In weak coupling theory, V is non-zero and 
assumed to be constant only within the narrow energy 
region |e — ei?| < lOc near to the Fermi energy ep, where 
LUc is the energy cut-off. 

Because of the pair-breaking effects of the strong 
exchange field in ferromagnets, we assume that only 
parallel-spin Cooper pairs can survive. Thus we set the 
p-wave antiparallel-spin gap function Aq = and retain 
the two gap functions with parallel-spin states ms = ±1, 
A-ti. The SC order parameter is assumed to have the 
following p-wave symmetry22, A±i(k) = (k^; + iky)A±. 

The Hamiltonian is treated via the Green function 
method within the mean-field theory framework. In ad- 
dition to the normal Green function ^^o.(k, r — r') = 
~{TTCk<j{T)cl^^{T')), the anomalous Green function de- 
scribing the pairing of electrons should be introduced, 
^cr(k, T — t') = {TrC\^a{T)c-\^a{T')) ■ Usiug the Standard 
equation of motion approach, the Green functions are 



derived to be 
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where the Pn are the Matsubara frequencies, and the FM 
and SC order parameters are respectively defined as 



^'^ = ^E((^k+)-(nk-)), 
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A±i (k) - - i ^ V^sc (k, k') jr± (k', r = 0). 



(3) 



All of the order parameters can be calculated using the 
above Green functions. They are found to satisfy 
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where = ek~fi±M, E±{k) - ^{ej^^f + |A±i(k)P, 
and f{E) is the Fermi function. The chemical poten- 
tial fi is determined from the equation for the number of 
electrons per unit volume, or particle density, 

^ 2i?_(k) 2E+{k) J' 

(6) 

which is equal to unity at half filling. 

Equations ^ , ^ and ([6]) with n — 1 comprise the self- 
consistent equations for the ferromagnetic superconduct- 
ing system. We solve the equations for the simple case of 
a spherical Fermi surface at half filling. It is convenient 
to solve these equations by converting the summations 
over k-space to continuum integrals over energy. 
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FIG. 1: The Stoner point Uc{V) as a function of the p-wave 
interaction strength V at T = 0. Inset: Enlargement of the 
region < F < 500. 

where ep^ — Jl ± M, e^-^ = e — ep_^, and E± ~ 
\J ie-^'^Y + si'^^ 6'|A±p. In the above equations, the unit 

t2 2/3 

of energy is rescaled by the factor . The di- 

mensionless interactions U and V are thus defined by 
jj ^ uiKf^Y^ and V = VlK^Y^, and the dimen- 
sionless energies 1f±, £, Wc, E±, A±, and Ji are defined 
analogously. The dimensionless temperature is defined 
by T = ksTi^^^^)-^. We choose ujc = O.Olei^, where 
EF is the dimensionless Fermi energy at M = T — 0. 

By solving the equations self-consistently, we can in- 
vestigate the interplay between the magnetism and the 
superconductivity in the coexisting state. This issue was 
discussed previously based on a similar framework, with 
the emphasis placed on the effects on the SC pairing due 
to the critical spin fluctuations in FM compounds'^^. The 
present work focuses on the reciprocal action, i.e., the in- 
fluence of the SC on the FM. 

According to Stoner theory, a Fermi gas can exhibit 
ferromagnetism only when the effective FM exchange is 
larger than the critical Stoner point. For a system de- 
scribed by Eq. (I), U represents the effective exchange 
interaction. In the absence of the p-wave SC interaction, 
V =_0, the dimensionless Stoner point Uc{0) « 12.76104. 
For V ^ 0, we calculate Uc{V). As shown in Fig. [U 
the T = Stoner point Uc{V) decreases as V increases, 
which implies that the p-wave Cooper pairing reduces the 
barrier to the onset of the magnetization of the Fermi gas. 
We note that V might be very small in a real system, so 
the enhancement effect of the superconductivity on the 
ferromagnetism may be very weak. The inset of Fig. [1] 
shows the details of Uc{V) in the region of small V, where 
the decreasing tendency of Uc{V) with increasing V still 
can be seen clearly. 

To further demonstrate the influence of the SC on the 
FM, we discuss the magnetization m = (n+) — (n_) as 
a function of F at T = 0. Here we use m ~ 2M/U 
instead of M to eliminate the dependence of U c upon 
V . As shown in Fig. [21 m{V) increases with V for each 




FIG. 2: Plots of the electronic magnetization m = {n+) — {n-) 
as a function of the p-wave interaction strength V at T = 
for fixed values of U. From larger to smaller m at fixed V, 
U = 12.8 (short dotted), 12.77 (dashed), 12.761 (solid), 12.743 
(dotted), 12.7 (dash-dotted) andl2.495 (short dashed). Inset: 
Enlargement of the region G <V < 500. 



given value of U. For U > Uc{0), m(0) is finite, since the 
system is spontaneously magnetized, and m{V) increases 
monotonically from m(0), eventually reaching unity at a 
finitc_F ^ 2300. For U <_C7c(0), however, m(V) = for 
V < Vc{U), and then rn(F) ^ increases sharply with V 
for V > Vc{U), eventually reaching unity at V > 2300. 
The critical value Vc{U) corresponds to the reduction in 
the Stoner point Uc{V) at which the onset of the ferro- 
magnetism is induced, as pictured in Fig. [TJ This is a 
second way in which the p-wave superconductivity can 
enhance the ferromagnetism. 

A similar effect was found in the ferromagnetic spin-1 
Bose gas which exhibits two phase transitions, the FM 
transition and Bose-Einstein condensation (BEC). The 
BEC temperature increases with FM couplings and, on 
the other hand, the FM transition is significantly en- 




FIG. 3: Plots of A+ (dashed) and A_ (solid) as functions of 
F at F = 12.77 and T = 0. Va is the value of V at which 
A_ has a maximum, and A_ a.t V — >~ 2300, the point 
at which m — > 1 in Fig. O 
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FIG. 4: Shown are plots of the order parameters M (dotted), 
A+ (dashed), and A_ (solid) as functions of T in the coex- 
istence state for V — 300. Af' (dash-dotted) is the magnetic 
order parameter when V = 0. (a) [/ = 12.79 > Uc{0) and 
> Tc+. (b)J7 = 12.7^> Uc{0) but < < T,+ . (c) 
U = 12.76 < Uc{0) but U > UciV). JThe ferromagnetism 
is induced due to the p-wave pairing (M 7^ 0) even though 
M' = 0. 



hanced due to the onset of the BEC^^i. Considering that 
triplet Cooper pairs behave somewhat like spin-1 bosons, 
a FM superconductor is analogous to a FM Bose gas. 

Figure [3] displays plots of the p-wave SC order param- 
eters, A± as functions of y at T = and U = 12.77, 
just above the V — Stoner point Uc{0). Although 
with increasing V, A_|_ rises monotonically, A_ initially 
rises, reaches a maximum at Va, and then decreases at 
an increasing rate until it vanishes discontinuously when 
m(V) = 1. For U = 12.77, miy) >_0 is sliown_by the 
dashed curve in Fig. ^ so that A_|_ > A_ for all V. Since 
m also grows with V, the mean number of spin-down elec- 
trons decreases with increasing V, vanishing when m — > 1 
at F « 2300, at and beyond which A_ 0. 

We now discuss the finite temperature properties of 
the system. We define M to be the magnetic order pa- 
rameter when V — 0, for which A± = 0. The T de- 
pendencies of the order parameters A±, M, and M are 
obtained numerically and shown for V = 300 and three 
different U cases in Fig. S) The order parameters be- 
come non-vanishing below their respective dimensionless 
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FIG. 5: Plots of the order parameters M (dotted), A+ 
(dashed), and A_ (solid) as functions of T in the coexistence 
state for F = 20 and t7 = 12.761 < Uc{0). 



transition temperatures Tc±, T^, and In each case, 
the SC order parameters A± increase monotonically with 
decreasing T below Tc±, respectively. In the FM super- 
conductor, Tc- < Tc+ and A_(T) < A+(T), as shown 
in Figs. 4(a), 4(b), and 4(c). In addition, m'(T) also 
increases monotonically with decreasing T for the ferro- 
magnet in the absence of any superconductivity, as de- 
picted in Figs. 4(a) and 4(b) for the respective cases 

17 > I7c(0) and t'^ > Tc+_and <t'^ <Tc+. How- 
ever, the T-dependence of M is non-trivial when p-wave 
superconductivity is present. In the first case pictured 
in Fig. 4(a), M(T) = m' {T) for t'^ > Tc+^ as in the 
absence of superconductivity. However, M (T) exhibits 
an upward kink at Tc+ below which A_|_ 7^ 0. Then, for 
Tc- < T < Tc+, M increases sharply with decreasing 
T, and exhibits a downward kink at Tc- below which 
A_ 7^ 0. Below Tc-, M{T) then decreases monotoni- 
cally with T . This case was discussed previously in a 
similar scenario^!. 

The case < Tc± not previously discussed is more 
interesting. Two examples of this case with V = 300 
are shown in Figs. Hl^b) and SJc). In Fig. 4(b), the 
magnetization M iov V = Q (and A± — 0) is so weak 

that < T^ < Tc-, but a non- vanishing V enhances 
the magnetization, M, causing the actual dimensionless 
Curie temperature to equal Tc+, below which both 
A+{T) and M{T) become discontinuously non- vanishing, 
signaling a first-order transition. Their behaviors for T < 
Tc+ — T„i are then qualitatively similar to that shown 
in Fig. 4(a), with A_(T) 7^ 0_for T < Tc-, cau_sing 
a downward kink in M{T) &iTc-, below which M(T) 
decreases monotonically with T. For the more extreme 
case when U < Uc{0) and t'„, = but 17 > Uc(V) 
depicted in Fig. 4(c), the behaviors of the three order 
parameters are very similar to that shown in Fig. 4(b). 

Considering that V is usually small in real systems, a 
case with V = 20 is checked, as shown in Fig. 5 where 
U is taken to be 12.761, slightly lower than Uc{0) but 
larger than C/c(20) « 12.7608. Fig. 5 looks very similar 
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to Fig. 4(c). _ 

Although we did not investigate the limit V 0+, 
the examples with V — 300 and F = 20 of the case 
T„j < Tc+ pictured in Figs. 4(b), 4(c) and Fig. 5 suggest 
that in FM superconductors, the actual Curie tempera- 
ture T„i is unlikely to ever be lower than the upper SC 
transition temperature Tc+, even if the FM order were 
extremely weak. In other words, these examples argue 
against the possibility of a FM T regime inside the p- 
wave triplet SC regime, with an actual Tm < Tc+. Anal- 
ogously, it was shown that the ferromagnetic transition 
never occurs below the Bose-Einstein condensation in the 
FM spin-1 Bose gaa^^. Moreover, the present results are 
to some extent consistent with the observed phase dia- 
grams of UGe2-^ and ZrZn2^, and with the theoretical 
discussion of Walker and Samokhin2^,^ who argued that 
the superconductivity only occurs within the FM region. 
In addition, this scenario is consistent with de Haas van 
Alphen experiments under pressure on UGe2^^. 

However, very recent experiments on UCoGe under 
pressure were interpreted as potentially having such a 
FM regime inside the SC regime near to the FM quan- 
tum critical poini*2^. However, the dc resistance and ac 
susceptibility measurements of T,„ and Tc+ could not de- 
termine if there were a FM region inside the SC one for 
pressures just below their extrapolated quantum critical 
pressure Pc, allowing for a first-order phase transition at 
the point when T,„ ~ Tc+, beyond which only a parallel- 
spin triplet state exists^*. Further experiments are en- 
couraged to determine if the FM and SC phase regimes 
with < T,„ < Tc+ at fixed pressure actually exist in 
UCoGe. 

As suggested by the results for the temperature depen- 
dencies of the order parameters, the FM superconduct- 
ing system shows multiple phase transitions, which can 
be determined experimentally from measurements of the 
specific heat. The specific heat at constant volume for 
our model can be calculated from 

C(T)=T^, 
^ ^ dT 

where the electronic contribution to the entropy S can 
be derived from 

^ = - E {/(^.)ln/(^a) 

k,o-=± 

+ f{E,)]\n[\- f{E,)]}. 

The specific heat was calculated previously based on 
a model of s-wave superconductivity coexisting with 
ferromagnetism^. For s-wave superconductors, there is 
only one SC transition temperature Tc, at which there is 
a jump in the specific heat at the second order transition. 
However, the case of a p-wave superconductor coexisting 
with ferromagnetism is more interesting. In Figs. 6(a) 
and 6(b), the results for the specific heat corresponding 
to the cases pictured in Figs. 4(a) and 4(b) for the order 
parameters are shown. For the case U > U c pictured in 
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FIG. 6: Plots of the electronic specific heat at constant volume 
as a function of T for (a) a case corresponding to Fig. Oa). 
The inset shows a transition from ferromagnetic to paramag- 
netic phase occurs at the Curie point Tm ~ 0.5. The dotted 
curve denotes the specific heat of the free electron gas; (b) a 
case corresponding to Fig. |4jb). The Curie point Tm = Tc+, 
at which the transition is first order. 



Figs. 4(a) and 6(a), there are three phase transitions at 
the temperatures Tc- < Tc+ < Tm- In Fig. 6(b), an ex- 
ample of the case < Tc+ when V — pictured in Fig. 
4(b) is shown. In this case with V — 300, there is a first- 
order phase transition at T„j — Tc+, and a second-order 
phase transition at Tc- ■ 

In conclusion, it is shown that p-wave triplet Cooper 
pairing can enhance the ferromagnetism in superconduct- 
ing ferromagnets. This enhancement is most prominent 
for the magnetic exchange interaction U very near to the 
Stoner point Uc{0), the critical value for the strength of 
the exchange interaction required for the onset of ferro- 
magnetism in the absence of the p-wave pairing interac- 
tion V. With finite V, Uc{V) is reduced and the fer- 
romagnetic order parameter increases in magnitude with 
increasing V. The temperature dependencies of the mag- 
netic and parallel-spin superconducting order parameters 
and of the specific heat are calculated. The results show 
that the Curie temperature is unlikely to ever be lower 
than the upper SC transition temperature, in agreement 
with pressure measurements on UGe2^- This feature 
also may be relevant to recent experiments on UCoGe^-. 
The temperature dependence of the specific heat exhibits 
two peaks for weak ferromagnetism in the coexistence 
state, with a first-order transition at the combined fer- 
romagnetic and upper p-wave SC transition, and a lower 
second-order p-wave SC transition. For strong ferromag- 
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netism, the specific heat exhibits three second-order tran- 
sitions. Our results support the possible coexistence of 
p-wa.ve superconductivity with a ferromagnetic state. 
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